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We study a free boundary problem modelling the growth of non-necrotic tumors with 
fluid-like tissues. The fluid velocity satisfies Stokes equations with a source determined 
by the proliferation rate of tumor cells which depends on the concentration of nutrients, 
subject to a boundary condition with stress tensor effected by surface tension. It is easy 
to prove that this problem has a unique radially symmetric stationary solution. By using 
a functional approach, we prove that there exists a threshold value 7* > for the surface 
tension coefficient 7, such that in the case 7 > 7* this radially symmetric stationary solution 
is asymptotically stable under small non-radial perturbations, whereas in the opposite case 
it is unstable. 
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1 Introduction 

In this paper we study the following free boundary problem modelling the growth of non-necrotic 
tumors with fluid-like tissues: 



Abstract 



Act 



f(a) in n(t) 



t > 





t > 




vAv + Vp 



_V(V-v) = 



in Q(t), t > 0, 





t > 




T(v,p)n = -7Kn on t > 0, 

V n = v • n on dQ(t), t > 0, 






v dx = 0, t > 




/ vxxdx = 0, t>0, (1.8) 
Jfl(t) 

n(o) = n , (i-9) 

where cr = cr(t, x), v = v(t,x) (= (vi(i, a;),«2(t, a;),«3(t, x))) and p = p(t,x) are unknown 
functions representing the concentration of nutrient, the velocity of the fluid and the internal 
pressure, respectively, / and g are given functions representing the nutrient consumption rate and 
tumor cell proliferation rate, respectively, which typically have the following forms respectively: 

f(a) = Xa, g(a) = »(a - a c ), (1.10) 

where A, (i and a c are positive constants, a c < a, and Q(t) is an a priori unknown bounded 
domain in M 3 representing the region occupied by the tumor at time t. Besides, v, a and 7 are 
positive constants, among which v is the viscosity coefficient of the fluid, 7 is the surface tension 
coefficient of the tumor surface, and a is the concentration of nutrient in tumor's host tissues, 
K, V n and n denote the mean curvature, the normal velocity and the unit outward normal, 
respectively, of the tumor surface dQ,(t), and T(v,p) represents the stress tensor, i.e., 

T(v,p) =f[V0v+(V®vf] _(p+Jlvv)I, (1.11) 

where I denotes the unit tensor. We note that the sign of the mean curvature k is defined such 
that it is nonnegative for convex hyper-surfaces. Without loss of generality, later on we assume 
that 

v = 1 and (7 = 1. 

Note that the general situation can be easily reduced into this special situation by using the 
rescaling a — > a/a, p — > p/v and 7 — ► j/u. 

Tumor growth modelling and analysis has attracted considerable attention during the past 
more than ten years. Most tumor models assume that the tumor tissue has the structure of a 
porous medium for which Darcy's law applies (see, e.g., [2], [3] and the references cited therein). 
For such tumor models, many interesting results of rigorous analysis have been obtained, for 
which we refer the interested reader to see [l]-[8], |18]-|22j. [26], [27] and the references cited 
therein. The tumor whose tissue does not have the structure of a porous medium but instead is 
more like a fluid was recently considered by Franks et al in the literatures [I2]-[T5], where some 
new models were proposed to mimic the early stages of the growth of ductal carcinoma in the 
breast. A basic feature of a ductal carcinoma in the breast in early stages is that it is confined 
to the duct of a mammary gland, which consists of epithelial cells, a meshwork of proteins, and 
extracellular fluid. In modelling, this leads to the replacement of the Darcy's law used in porous 
medium structured tumor models by the Stokes equations. See p2]-[15] for details. The models 
of Franks et al |12j-[15] have been concisely reformulated by Friedman in [16] (see also [17]). 

The problem (1.1)— (1.9) above is a simplification of the tumor model proposed in [16] . The 
simplifications are made in two aspects. First, the model in [16] contains a system of nonlinear 
hyperbolic conservation laws with source terms describing the movements and interchanges of 
three different species of cells: proliferating cells, quiescent cells and necrotic cells. In this paper 



we only consider one species of proliferating cells, so that no hyperbolic conservation laws appear 
in (1.1)-(1.9). Second, in [16] the equation for a is of the following evolutionary type: 

d t a = Act - f(a) in fi(t), t > 0, 

where / is as given in (1.10), but in this paper the stationary form (1.1) is considered. All these 
simplifications are made for the purpose to make the analysis simpler. If either one of the above 
two aspects of simplifications are not made, then the model will be much more complicated to 
analyze, and new mathematical techniques have to be employed. We leave it for future work. 

In [16] Friedman established local wellposedness in Holder spaces of his model. Meanwhile, 
he proved that in the special case that the tumor contains only one species of cells (i.e., the tumor 
contains only proliferating cells), there exists a unique radially symmetric stationary solution. 
Based on these results, a number of interesting questions are raised in [16] (see also [17]), one of 
which is as follows: Is this radially symmetric stationary solution asymptotically stable under 
non-radial perturbations? A heuristic result toward an answer to this question was obtained by 
Friedman and Hu in [19], where they proved that this radially symmetric stationary solution is 
linearly asymptotically stable for small /x/7, i.e., there exists a threshold value (/i/7)* such that 
if we denote by (cr s ,v s ,p s , £l s ) this stationary solution, then in the case /i/7 < (/i/7)* the trivial 
solution of the linearization at (a s , v s ,p s , f2 s ) of the original problem is asymptotically stable. 
Moreover, they also proved that in the case /i/7 > (/x/7)* the radially symmetric stationary 
solution is unstable. We also refer the interested reader to see |18] for the study of existence of 
non-radial stationary solutions. 

The purpose of this paper is to prove that, at least for the simplified model (1.1)— (1.9), the 
answer to the above question is yes for large 7 but no for small 7. To this end, we shall use a 
functional approach inherited from the references [7], [10], [26] and [27] to study the problem 
(1.1)-(1.9), namely, we shall first reduce the problem (1.1)— (1.9) into an evolution equation 
containing merely the function p describing the free boundary <9f2(i), which can be considered 
as a differential equation in certain Banach space. We shall prove that this differential equation 
is of the parabolic type. Next we use the geometric theory for parabolic differential equations 
in Banach spaces (see [1] and [2J]) to study the stability of the stationary solution. Since 
our discussion does not depend on the specific linear forms of the equations (1.1) and (1.2), 
throughout this paper we shall not consider the specific forms of / and g given by (1.10), but 
instead assume that they are general smooth functions satisfying the following assumptions: 

(Al) f G C°°[0, 00), f'(a) > for a > and /(0) = 0. 

(A2) g G C°°[0, 00), g'(a) > for a > and g(a c ) = for some a c > 0, 

(A3) a c < 1. 

To give a precise statement of our main result, let us first introduce some notations. 

Given a bounded domain !!C1 3 and two numbers m G N and 6 G (0, 1), we denote by 
h m+ (Q) the so-called little Holder space on O of index m + 6, which is, by definition, the closure 
of C°°(Q) in the usual Holder space C m+e (Q). Similarly, given a smooth hypersurface T in M 3 , 
we denote by h m+e (T) the closure of C°°(T) in C m+0 (T). 

It can be easily shown (see Theorem A in Appendix A) that under Assumptions (Al)— 
(A3), the problem (1.1)-(1.8) has a unique radially symmetric stationary solution. Later on 



we use the same notation (<r s , v s ,p s , S ) as before to denote this radially symmetric stationary 
solution of the problem (1.1)-(1.8). Note that this means that there exists R s > such that 
Q s = { x G M 3 : \x\ < R s } and 

X 

a s (x)=a s (r), v s (x) = v 8 (r)-, p s (x) = p s (r) for x G O s , (1-12) 

r 

where r = \x\ and t> s represents a scalar function. Clearly, a coordinate translation of a solution 
of (1.1)-(1.8) is still a solution of it. Thus, for any xq G M 3 , we denote by (owii v b ]>P[cco]' ^teol) 
the stationary solution obtained by the coordinate translation x — > x + xo of the stationary 
solution (a s ,-v s ,p s ,Q, s ). Given p G C 1 (c?r2 s ) with HpHc 1 ^^) sufficiently small, we denote by Cl p 
the domain enclosed by the hypersurface r = R s + p(co), where u> G 00^. It is obvious that for 
xq G M 3 such that \xq\ is sufficiently small, there exists a smooth function pr^i on d£l s such 
that 0[ So ] = Slp^ y Since we shall only be concerned with small perturbations of the stationary 
solution (o~ s , v s , p s , Q s ), it is natural to assume that the domains £l(t) and S7o in (1.1)— (1.9) are 
small perturbations of £l s . It follows that there exist functions p(t) (= p(u,t)) and po (= po(k>)) 
on d£l s such that 0(t) = Fl p (t) and f^o = Q po - Using these notations, the initial condition (1.9) 
can be rewritten as follows: 

p(u, 0) = Pq(w) for u> G d^l s . (1-13) 

The solution (a,v,p,£l) of the problem (l.l)-(l. 9) will be correspondingly rewritten as (a,v,p, p), 
and the radially symmetric stationary solution (a s ,v s ,p s ,£l s ) will be re-denoted as (er s , v s , p s , 0) . 

The main result of this paper is the following theorem: 

Theorem 1.1 Assume that Assumptions (Al)-(A3) hold. For given m G N, m > 3, and 
< 6 < 1, we /iaue t/ie following assertion: There exists a positive threshold value 7* such that 
for any 7 > 7*, i/ie radially symmetric stationary solution (a s ,v s ,p s ,0) is asymptotically stable 
in the following sense: There exists constant e > such that for any po G h m+d (dQ s ) satisfying 
\\po\\c m + e (dn s ) < e > problem (1.1)— (1.9) /ias a unique solution (a,v,p,p) for all t > 0, and 
there are positive constants to, K independent of the initial data and a point xq G M 3 uniquely 
determined by the initial data, such that the following holds for all t > 0: 

Ik(')t) - a [x ]\\c m + e (n(t)) + ll v (->*) ~ v [x ]\\c m - 1 + e (Q(t)) 
+||p(-. *) - P[a> ] lie— 2 + e (n(t)) + ll/ 5 !-' *) - P[* ] llc™+ e (a^) < Ke-"*. (1.14) 

For 7 < 7* f/ie stationary solution (a s ,v s ,p s ,0) is unstable. □ 

It is interesting to compare this result with the corresponding result for the porous medium 
structured tumor model obtained by Cui and Escher in [8], where it is proved that, for the 
porous medium structured tumor model, there exists a threshold value for the surface tension 
coefficient 7, which we denote as 7*, such that the unique radially symmetric stationary solution 
is asymptotically stable if 7 > 7*, but unstable if 7 < 7*. We shall show that 7* > 7*. This 
implies that radially symmetric stationary solution is more stable for a tumor whose tissue has 
a porous medium structure than a tumor whose tissue is more like a fluid. See Lemma 3.5 for 
the proof of the assertion that 7* > 7*. 

The structure of the rest part is as follows. In Section 2 we first convert the problem into 
an equivalent initial-boundary value problem on a fixed domain by using the so-called Hanzawa 



transformation, and next further reduce it into a scalar equation containing the single function 
p, which can be regarded as a differential equation in the Banach space /i m ~ 1+6l (§ 2 ). We shall 
also prove that this equation is of the parabolic type. In Section 3 we study the linearization 
of (1.1)— (1.8) at the radially symmetric stationary solution, and study the spectrum of the 
linearized operator. In the last section we give the proof of Theorem 1.1. 



2 Reduction of the problem 

In this section we reduce the problem (1.1)-(1.9) into a differential equation in a Banach space. 
For simplicity of the notation, later on we always assume that R s = 1. Note that this assumption 
is reasonable because the case R s ^ 1 can be easily reduced into this case after a rescaling. It 
follows that 

!l s = B 3 = {i;GK 3 : |x| < 1} and dVL s = § 2 = {x G E 3 : |x| = 1}. 

Let m and 9 be as in Theorem 1.1. For u G h m+e (M 3 ), we denote by tr§2(u) the trace of 
u on § 2 , i.e., tr§2(u) = u|§2. We know that tr§2(u) G h m+e (§> 2 ) and the operator tr S 2 : u — > 
tr§2(u) from h m+0 (M ) to h m+ (S 2 ) is linear, bounded and surjective. Let II be a bounded 
right inverse of it, i.e., II G L(h m+e (S 2 ), h m+e (M 3 )) and tr§ 2 (n(u))) = u for any u G h m+e (S 2 ). 
Let E G L(C m+e (M 3 ),BUC m+e (R 3 )) be an extension operator, i.e., E has the property that 
E(u)(x) = u{x) for any u G C m+e (E 3 ) and x G B' 3 . Here BUC m+e (R 3 ) denotes the space 
of all C rn functions u on M 3 such that u itself and all its partial derivatives of order< m are 
bounded and uniformly 6-th order Holder continuous in M 3 . We denote Hi = E o H. Then 
clearly n x G L(h m+e (§ 2 ), h m+e (R 3 )), where h m+e (R 3 ) represents the closure of BUC°°(R 3 ) in 
BUC m+9 (R 3 ). Hence there exists a constant C > such that 

\\^i(p)\\BUCm+e m < CoWpWcm+e^ for p G /^(S 2 ). (2.1) 

Take a constant < 5 < min{l/6, l/(3Co)} and fix it, where Co is the constant in (2.1). We 
choose a cut-off function x £ C°°[0, oo) such that 

We denote 

Given p G 0™ +6> (§ 2 ), we define the Hanzawa transformation <E> p : R 3 —>■ R 3 as follows: 

<S> p (x) =x + x(r-l)Tli(p)(x)^ for iGR 3 . (2.3) 

Using (2.1) and (2.2) we can easily verify that & p is a h m+e diffeomorphism from R 3 onto itself, 
i.e., <E>p G Diff m+6, (1R 3 ,1R 3 ), and each component of $ p and belongs to h m+0 (R 3 ). Later on 
we write <3? p G Diff^ +e (R 3 , IR 3 ) to indicate this fact. We define <j) p = ® P \ s2 , and denote 

n p = $ p (b 3 ), r p = dn p = $ p (s 2 ) = im(^). 



Clearly, 

4>p(ui) = [1 + p{u))]u) for uj G 8 2 . 

This implies that x G T p if and only if there exists w £ § 2 such that x = [1 + p(u;)]u;. Thus, 
in the polar coordinates (r,u>) of M 3 , where r = \x\ and w = x/\x\, the hypersurface T p has the 
following equation: r = 1 + 

Next, given p G C([0, T], C^ t+9 (§ 2 )), for each t G [0, T] we denote 

r p0) = r P (t) and fy>(*) = tt P {ty 

Since our purpose is to study asymptotical stability of the radially symmetric stationary solution, 
later on we always assume that the initial domain Q.q is contained in a small neighborhood of 
$l s = B 3 . It follows that there exists p G Of +e (S 2 ) such that T = dfl = F po . 

Let p be as above, and let $ p be the i-th component of <J> P , i = 1,2, 3. We denote 

[D* p ] ir .= di& p = -^, a^ j (x) = [D^ p (x)]- j 1 (i,j = 1,2,3), 

G p (x) = det(D<S> p (x)) for x G M 3 , 

= \</> P \ 2 yJl + \V u <i> P \ 2 for u, G S 2 , 

where V w represents the orthogonal projection of the gradient onto the tangent space 
TL^S 2 ) 1 ). Here and hereafter, for a matrix A we use the notation Aij to denote the element of 
A in the (i,j)-th position. Since $ p G Diff^ +e (M 3 , M 3 ), we have a?- G /i m - 1+f? (lR 3 ), i,j = 1,2,3, 
G p G h m ~ 1+9 (M. 3 ), and G h m ~ 1+e (S 2 ). We now introduce four partial differential operators 
A{p), B(p), B(p)- and B{p)® on E 3 as follows: 

A{p)u{x) = a^-{x)dj{a p ik {x)dku{x)^ for scalar function u, (2.4) 

B(p)u(x) = ( y a1j(x)dju(x),a2j(x)dju(x),a^j(x)dju(x)^ for scalar function u, (2.5) 
B(p) ■ v(x) = a1-{x)djVi{x) for vector function v = (^1,^2,^3)- (2-6) 
S(p) (8) v(x) = (a^ k (x)dkVj(x)) for vector function v = (fi, f 2, f 3). (2.7) 

Here and hereafter we use the convention that repeated indices represent summations with 
respect to these indices, and dj = d/dxj, j = 1, 2, 3. Obviously, 

A(p) G L(h m+e (B 3 ), h m ~ 2+e (B 3 )), B{p) G L(h m+e {M 3 ), (h m ~ 1+e {M 3 )) 3 ), 

B(p)- G L{(h m+e (E 3 )) 3 ,h m - 1+e {M 3 )), B{p)® G L((h m+9 (M 3 )) 3 , (h m - 1+e (M 3 )) 3x3 ) . 
The definitions (2.4)-(2.7) can be respectively briefly rewritten as follows: 

A{p)u = (A(u o o $ p , B{p)u = (V(u o ^)) o <D p , 



1 - ) In the coordinate w = = (sini9 cos sin ■& sin <p, cos #) (0<$<7r, 0<<^< 27r) of the sphere we have 

V w f(u>) = (cos 1? cos cos 1? sin — sin#)9#/(u;(#, <^)) H — ^"~^(~ sin^, cos 99, 0)d v f(Lj('d, ip)). 

df 1 
Note that V x f = tt-w + -V w /. 

or r 



B(p) • v = (V • (v o S- 1 )) o $ p , ® v = (V ® (v o S- 1 )) o $ p . 

As in |16j we introduce the following vector functions: 

wi(z) = (0, x 3 , -x 2 ), w 2 (x) = (-x 3 , 0, xi), w 3 (x) = (x 2 , -xi, 0). 

Then clearly v x x = (v • W],v • w 2 , v • W3). 

Let n and k be respectively the unit outward normal and the mean curvature of T p (see 
(1.5)). We denote 

n p (x) = n(0p(x)) and K p (x) = k(4> p (x)), for x G S 2 . 
A direct computation shows that 



n„ x 



x- [(ZJ^(x))- 1 ] 
\x-[(DS> p (x))-i] 



a p ij (x)x j e l 



|af-(x)x i e i 



= (Si(20,H?(a?),Hp(:c)), 



(2.8) 



where 



and 



ei = (1,0,0), ea = (0,1,0), ea = (0,0,1), 
«p(s) = ^'(^HH^x). 



(2.9) 



We remind the reader to notice that, obviously, [p -» $ p ] G C°°(0^ +9 (S 2 ), Diff™ +9 (1R 3 , 
Thus we have 

*,i = 1,2,3, 



V 


- <•] 


G C°°(0™ +e 


(S 2 ),/i m - 1+ 


9 (M 3 )), 


[p- 




G C°°(0™ +6 


(S 2 ),/i m - 1+ 


e (K 3 )), 


[p- 




G C°°(0™+ £ 


(§ 2 ),/i m ~ 1+ 


*(§ 2 )), 


[p- 


-> « p ] 


G C°°(0^ +e 


(S 2 ),/i m " 2+£ 


(§ 2 )), 


>- 




G C°°(0™ +e 


(S 2 ),(/i m ~ H 


- e (S 2 )) 3 



(2.10) 



Finally, for a, v and p as in (1.1)-(1.8), we denote 

a = a o <3? p , v = v o p = po $ p . 



We also denote w 



,1 



w,-o$„, j = 1,2,3. 



Using these notations, we see easily that the Hanzawa transformation transforms the equa- 
tions (1.1)— (1.5), (1.7) and (1.8) into the following equations, respectively: 



A(p)a = /(or) in B 3 , t > 0, 
B(p) ■ v = 5(0?) in B 3 , t > 0, 
-A(p)v + B(p)p - l -B{p){B{p) • v) = in B 3 , 



t > 0, 



(7 = 1 on S 2 



t > 0, 



T p (v,j?)n p = -7K p n p on S 2 , f > 0, 



(2.11) 

(2.12) 

(2.13) 

(2.14) 
(2.15) 



/ v(x)G p (x)dx = 0, t > 0, (2.16) 

J\x\<l 

/ v(x) • w p Jx)G p (x)dx = 0, j = 1,2,3, t > 0. (2.17) 

^|x|<l 

Here T p (v,p) = [B(p) <g> v + <g> v) T ] - [p + (2/3)B(p) ■ v]I. 

In what follows we rewrite (1.6) into an explicit equation expressed with the function p = 
p(uj,t). Let ip p (x,t) = r — 1 — p(uj,t), where r = \x\ and to = x/\x\. Then x G T p (i) if and only 
if ip p {x,t) = 0. It follows that the normal velocity of T p (t) is as follows (see |10|): 

Fn(a?,t)= J^Am for x G r p (t), t>0. 

Moreover, n(x,i) = V x ip p (x, t)/\V x ijjp(x, t)\. Hence (1.6) can be rewritten as follows: 

d t p(u),t) = v(x,t) ■ V x ip p (x,t) for x G T p (t), t > 0, 
dip 

where uj = x/\x\. Since V^Vv = —-^-u)-\-r~V UJ ibp = w-r^V^p, we see that after the Hanzawa 

or 

transformation, this equation has the following form: 



V w p(a»,t) 

• \uj- 

Finally, we rewrite (1.13) as follows: 



d tP (u, t) = v(w, t)-[u- - ' ^ ] for we§ 2 , t > 0. (2.18) 



p(u,0) = Po{u) for w 6 § 2 . (2.19) 
In summary, we have the following preliminary result: 

Lemma 2.1 If (a,\,p, p) is a solution of the problem (1.1) -(1.9), then by letting a = a o& p , 
v = v o <£p andp = poQ p , we have that (<x,v, p, p) is a solution of the problem (2. 11) -(2. 19). 
Conversely, If (a, v,p, p) is a solution of the problem (2.11)-(2.19) ; i/ien fry letting a = <jo<3?~ 1 , 
v = vo^ 1 andp = po&~ 1 , we have that (a,v,p, p) is a solution of the problem (1.1) -(1.9). □ 



In the sequel we further reduce the problem (2.11)-(2.19) into a scalar equation containing 
the single unknown p. The idea is to first solve the system of equations (2.11)-(2.17) to get a, 
v and p as functionals of p, and next substitute v obtained in this way into the equation (2.18). 

We first consider (2.11) and (2.14). We have: 

Lemma 2.2 Let 5 be sufficiently small. Then, given p G 0™ +e (S 2 ), the boundary value 
problem 

A(p)a = /(of) in B 3 , a = l on § 2 (2.20) 
has a unique solution a = TZ(p) G h m+e (M> 3 ) which satisfies < a < 1. Moreover, we have 

n £ C°°(0™ +e (S 2 ), h m+e (B 3 )). (2.21) 



Proof: See Lemma 3.1 of [7] ■ D 

Next, for given p G 0™ +6l (§ 2 ) we consider the following boundary value problem: 

B(p) -v = <p in B 3 , (2.22) 

-A(p)v + B(p)p = g in B 3 , (2.23) 
T„(v,p)£„ = h on § 2 , (2.24) 

v(x)G p (x)dx = 0, (2.25) 

x|<l 

f v(x)-w p Jx)G p (x)dx = 0, j = l,2,S, (2.26) 
./|a;|<l 

where tp G /i™"*- 14 ^! 3 ), g G (/i m ^ 2 ^(B 3 )) 3 and h G (/^-^(S 2 )) 3 for some 0< fc <m-2. 

Lemma 2.3 Lei 5 be sufficiently small and let p G 0™ +e (S 2 ) be given. A necessary and 
sufficient condition for (2.22) -(2.26) to have a solution is that ip, g and h satisfy the following 
relations: 

[ {z{x)-\B{p) V {x))-^{x)G p (x)dx+ [ h(x)-^(x)H p (x)dS x = 0, j = 1,2,3, (2.27) 

J\x\<l J\x\=l 

[ (g(x) - ^(pMz)) • ejGpWdx + / h(x) • ejH p (x)dS x = 0, j = 1, 2, 3. (2.28) 

J\x\<l 3 J\x\=l 

If this condition is satisfied, then (2.12)-(2.26) has a unique solution (v,p) G (/i m_fe+£ '(B 3 )) 3 x 
/i m ' fc - 1+e (B 3 ). 

Proof: Integrating by parts and employing the divergence theorem we see that for any v, 
w G (C 2 (f2 p )) 3 and p G C 1 (r2 p )) there holds the following integral identity (cf. [TT]): 



- / [V®v|(V« v) T ] • [V w + (V w) T ]dx -A [p + -(V ■ v)]V • wcfe 
= / [- Av + Vp V(V • v)] • wdx + T(v, p)n • wdS : 



p 



(2.29) 



Here, for two matrix A and B we use the notation A ■ B to denote the inner product of A 
and -B, i.e., A ■ B = AijBij. By making the Hanzawa transformation, we see that for any v, 
w G (C 2 (B 3 )) 3 and p G C 1 ^ 3 ) there holds 

\ [ [B(p)^v+(B(p)^) T ]iB(p)^w+(B(p)®^) T ]G p dx- ( [p+\B{p)-w]B{p)-^G p dx 

^J\x\<X J\x\<X 

= ! [~A{p)v + B(p)p- \b{ P ){B{ P ) ■ v)] • wG p dx + / T p (v,p)n p • wff,^. 

^|s|<l J|x|=l 

(2.30) 

Besides, clearly V (g> + (V ® w j) T = 0, V • Wj = 0, j = 1, 2, 3, which yield, after the Hanzawa 
transformation, that 



B(p) ® + 8) wp T = 0, B(p) • = 0, j = 1, 2, 3. 



Hence, if (v,p) is a solution of (2.22)-(2.26), then by takeing w = (j = 1,2,3) in (2.30), we 
see that (2.27) holds. Similarly we have (2.28). This proves the necessity of (2.27) and (2.28). 

Next we assume that the conditions (2.27) and (2.28) are satisfied, and proceed to prove 
that there exists a unique solution to the problem (2.22)— (2.26). We first prove uniqueness of 
the solution. Let (vi,^) and (v2,P2) be two solutions of (2.22)-(2.26). Then v = vi — V2 and 
p = Pl —p2 satisfy the corresponding homogeneous equations. Thus, by letting w = v in (2.30), 
we get 

/ 0(j>) ® v + (B(p) ® v) T ] • [B(p) ®w + (B(p) ® v) T ]G p (x)dx = 0, 

J\x\<l 

so that 

B(p) ® v + (B(p) ® v) r = in B 3 . 
This combined with (2.25) and (2.26) yields, by the Korn inequality 



U| l(i?!(^)) 3 - C ' 1 ll S ( U )ll(L 2 (n))3x3 
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/ u x xdx 
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(2.31) 



where S(u) = V ® u + (V tg> u) T (cf. Proposition 8.1 of [25]), that v = 0. From this it follows 
immediately that also p = 0. Hence the solution is unique if it exists. 

To prove existence we denote, for a given < k < m — 2, 

-fc-l+6»/W 3 \ ttd6 



X = (/l r ' 



- fc+e (B 3 )) 3 X // 
-fc-2+6»/Tffi 3 ^3 



Y = h m - k - 1+e (M 3 ) x (/i m - fc - 2+e (B a )) 3 x (h m - k - 1+e (S 2 )) 3 x E 3 x M 3 , 

and regard 0™ +6l (§ 2 ) as an open subset of the Banach space h m+e (§ 2 ). For every p E 0™ +6> (S 2 ) 
we define a linear operator £(p) : X — > Y as follows: 



£(p)C7 



• v 

-■4(p)v + tf(p)p + 
T p (v,p)S p 

v(x)G p (x)dx 

x|<i 

v(x) • w p Ax)G p {x)dx 

\x\<\ 



for U = (y,p,C) £X, 



where L is the linear operator from R° to (/t' 
6 3 wg for C = (a, 6 X , 6 2 , &a) E M 3 x R x R 



m-k-2+6 



)) 3 defined by /(C) = a + bi w? + 6 2 w^ + 



By (2.10) it is clear that 



£ E C°°(0™ +e (S 2 ),L(X,Y)), 



and 



C(0)U = 



where Iq is the linear operator from R 6 to (/i m ~ fe ~ 2+6l (B 3 )) 3 defined by = a + 61 wi + 62W2 + 
63W3 for ( = (a, 61,62,63) G R 3 x R x R x R. From the proof of Theorem 2.1 of [16] we see 
that C(0) is an isomorphism from X to Y (cf. also Lemma A.l of [9j). Since all isomorphisms 
from X to Y forms an open set in L(X,Y), we conclude that for 5 sufficiently small, C(p) 
is also an isomorphism from X to Y for any p G 0™ +e (S 2 ). This particularly implies that 
given ip G /i m - fe - 1+e (B 3 ), g e (/i m - fc - 2+e (B~ 3 )) 3 and h G (/i m - fe - 1+e (B 3 )) 3 , there exist unique 
v G (/i m " fc+e (B 3 )) 3 , p G /i m - fe - 1+9 (B 3 ) and ( 6 R 6 such that they satisfy (2.22), (2.24)-(2.26) 
and 

-A(p)v + + l p (() = g in B 3 . (2.32) 

We claim that ( = 0. Indeed, taking w = w^e, in (2.30) and using (2.22), (2.24), (2.27), (2.28) 
and (2.32), we get 

( l p (C)-^G p dx = 0, ( l p (()-e 3 G p dx = 0, j = 1,2,3. 

i|x|<l 7|x|<l 

From these relations we can easily show that if p = then £ = 0. By continuity (a small 
perturbation of a nonsingular matrix is still nonsingular) , this implies that if 5 is sufficiently 
small then for any p G 0™ +6 *(§ 2 ) we also have £ = 0. Hence our claim is true. It follows that 
(v,p) is a solution of (2.22)-(2.26). This completes the proof of Lemma 2.3. □ 

Lemma 2.4 For the solution of the problem (2.22) -(2.26), we have v = V(p)ip + Q(p)g + 
R(p)h, where 

m—2 

re p c™(o™ + \s 2 ),L(h m - k - 1+e (M 3 ),(h m - k+e (n 3 )) s ), 

k=0 
m-2 

< QG p| C°°(0™ +e (S 2 ),L((/i m ^ 2+e (B 3 )) 3 ,(/i m - fc+e (B 3 )) 3 )), 

fe=0 
m-2 

RG P C oo (O™ +0 (S 2 ),L((/i m - fc - 1+e (S 2 )) 3 ,(/i m - fc + e (B 3 )) 3 )). 

fc=0 

Proof: Let notations be as in the proof of Lemma 2.3. We denote by I\, I2 and I3 
the natural embedding operators from h m ~ k - 1+e (if), {h m - k ~ 2+e '(B 3 )) 3 and (§ 2 )) 3 , 

respectively, into Y, and by J7" the projection operator from X onto (/i m-fc+9 (B )) 3 . Then by 
letting 

V(p) = JoC(p)- 1 ol 1 , Q(p) = JoC(p)- 1 ol 2 , K(p)=JoC(p)- 1 ol 3 , 



V -v 

-Av + Vp + / (C) 
T(v,p)n 

v(x)da; 

x|<l 
v(x) x xcix 

|as|<l 



for 17= (v,p,C) G X, 



we immediately see that the desired assertion follows. □ 

The system of equations (2.12), (2.13) and (2.15)-(2.17) can be rewritten in the form of 
(2.22)-(2.26), with 

V = g(&), g = -^B{p)g(a), h = -jK p n p . (2.33) 

We assert that (2.27) and (2.28) are satisfied by these functions. Indeed, since g = ^B(p)ip, this 
assertion follows if we show that 

f h(x)-w P (x)H p (x)dS x = 0, [ h(x)- ej H p (x)dS x = 0, j = 1,2,3. (2.34) 

J\x\=l J\x\=l 

Let A-p p be the Laplace-Beltrami operator on T p . Then we have n(x)n(x) = —Ay p x for x G T p 
(cf. [9], [25]). Since Ar p is a symmetric operator in L 2 (T p ,dS x ), we see that 

/ (Ar p xi ■ xj - A Tp Xj ■ Xi)dS x = 0, i,j = 1,2,3. (2.35) 
Jv 



Thus 



/ K p n p ■ w P H p dS x = ™ • WjdS x = - A T x • wjdS x = 0, j = 1, 2, 3. 
J\x\=i Jr a Jr a 



Similarly we have 

/ K p n p ■ ejHpdS x = kii- ejdS x = - A Tp x ■ ejdS x = 0, j = 1, 2, 3. 
J\x\=i Jr p Jr p 

This verifies (2.34). 

Now, given p £ 0™ +e (§ 2 ), we first use Lemma 2.2 to solve the equations (2.11) and (2.14). 
This gives a = K(p) G h m+e (B 3 ). Next we use Lemma 2.3 to solve equations (2.12), (2.13) and 
(2.15)-(2.17). Note that with tp, g and h given in (2.33), we have <p 6 h m+9 (B 3 ) C h m - 2+d (\f), 
g G (h m ~ 1+e (B 3 )) 3 C (/i m ~ 3+e (B 3 )) 3 , and, by (2.10), h G (h m ~ 2+e (S 2 )) 3 . Hence, by Lemma 2.3 
(with k = 1) it follows that these equations have a unique solution (v,p) G (/i m_1+e '(B 3 )) 3 x 
/i m ~ 2 + e (B ). Moreover, since a = 1Z(p), by Lemma 2.4 we have 

v = V{ P )g{n{p)) + l -Q{p)B{p)g(K(p)) - 7 R(p)(lC(p)tf(p)). (2.36) 

where fC(p) = k p and M{p) = n p . We note that 

K. G C°°(Of +9 (S 2 ), /i m ~ 2+e (S 2 )), AT g C°°(0™ +9 (S 2 ), (h m - 1+e {§ 2 )) 3 ). (2.37) 

Substituting the expression of v in (2.36) into (2.18), and introducing the operator Q : 0™ +6l (§ 2 ) - 
h m - 1+0 (§ 2 ) by 

Q(p) = tr §2 [V(p)g(K(p)) + ±Q(p)B(p)g(K(p)) - 7 R(p)(lC(p)tf(p))} ■ [u - (2.38) 

(for p G 0™ +e (§ 2 )), where as before u represents the variable in § 2 , we see that the problem 
(2.11)-(2.19) is reduced into the following initial value problem for a differential equation in the 
Banach space h m ~ 1+e (§> 2 ): 

r*.-ew, *>o, (239) 

I p\t=o = po- 



Here Q is regarded as a unbounded operator in h m 1+6 (S 2 ) with domain 0™ +6, (§ 2 ). 
We summarize: 

Lemma 2.5 Let (a,v,p,p) be a solution of the problem (2.11) -(2.19). Then p is a solution 
of the initial value problem (2.39). Conversely, if p is a solution of the initial value problem 
(2.39), then by letting a = IZ(p) and (v,p) = P£(p)~ 1 (g(a), ^B(p)g(a), — 7K p n p , 0, 0), where 
P denotes the projection from X = (/i m - 1+9 (B 3 )) 3 x h m - 2+9 (M 3 ) x R 6 onto (h m - 1+e (B 3 )) 3 x 
h m ~ 2+e (M 3 ), we have that (a,v,p,p) is a solution of (2.11) -(2.19). □ 

From (2.10), (2.21), (2.37) and Lemma 2.4 we see that 

Q G C°°(0^ +e (§ 2 ), h m - 1+e {§ 2 )). (2.40) 

In the sequel we prove that if 5 is sufficiently small then for any p G 0™ +e (S 2 ), DQ(p) is a 
infinitesimal generator of an analytic semigroup in /i m_1+e '(S 2 ) with domain h m+9 (S 2 ), so that 
the differential equation in (2.39) is of parabolic type. Here and in what follows, the notation 
D- represents Frechet derivatives of smooth operators from h m+e (S 2 ) to /i m ~ 1+e (§ 2 ). 

We first note that the mean curvature operator KL has the following expression 

K(p)=}C 1 ( P )p + Ko(p), (2.41) 

where, for each p, )Ci(p) is a second-order linear elliptic partial differential operator on §> 2 with 
coefficients being functions of p and its first-order derivatives, and /Co is a first-order nonlinear 
partial differential operator on S 2 , so that 

Ki G ^(O^+^S 2 ), L(h m+e {S 2 ), h m - 2+e {S 2 ))), /C G C°°(0™ +9 (§ 2 ), h m - 1+e {S 2 )). (2.42) 

(see [7] and [TP]). Substituting (2.41) into (2.38) we see that 

Q{p) = Q2(p)[p,p] + Qi(p)p + Qo(p), (2.43) 

where, for each p, Q 2 {p) is a bilinear operator, Qi(p) is a linear operator, and Qq is a nonlinear 
operator; they are respectively defined as follows: 

Q2( P )[m,m] =jtr$2{R(p)[K: 1 (p) m rt(p)]}- 

Q^r, =-t^{V(p)g(TZ(p)) + l -Ct{p)B{p)g{lZ{p)) - 7 R(p) [JC (p)tf(p)] } " ^ 

-7trs2{R(^[X:i(p)^(p))]}-a;, 
Q (p) = tr §2 {V(p)g(n(p)) + l -Ct(p)B(p)g{lZ{p)) - 7 R(p) [K* {p)ti{p)\ } ■ u. 

We note that, by (2.21), (2.37), (2.42) and Lemma 2.4 we have 

Q 2 G C^(0™ +e (S 2 ),£L(/i m+e (S 2 ) x h m+e {§ 2 ),h m ~ 1+e {§ 2 ))), 
Qx G C°°(0^ +e (S 2 ),L(/i m + e (S 2 ),/i m - 1 + e (S 2 ))), 
Qo G C° o (0™ +e (S 2 ),/ i m+e (S 2 )), 



where BL{- x •, •) denotes the Banach space of all bilinear operators with respect the correspond- 
ing spaces. 

Given two Banach spaces Eq and E\ such that E\ is continuously and densely embedded 
into Eq, we denote by H(Ei, Eq) the subset of all linear operators A G L(E%, Eq) such that —A 
generates a strongly continuous analytic semigroup on Eq. 

Lemma 2.6 -DQ(0) G H{h m+9 (§ 2 ), h m - 1+0 (S 2 )). 

Proof: For any p G Of +e {S 2 ) and 77 G h m+e {S 2 ) we have 

DQ(p) V = Q 2 (p)[r],p] + Q 2 (p)[p,r?] + [DQ 2 (p) V ][p, p] + Qi(p)r? + [DQx(p)r)]p + DQ (p)r]. 

In particular, 

DQ(0)?7 = 21(0)7/ + 1)20(0)7/ for 77 G /i m+e (§ 2 ), (2.44) 

i.e., 1)2(0) = Qi(0) + DQ (0). We note that Qi(0) G L(/i m+e (§ 2 ), /i m - 1+e (§ 2 )) and DQ (0) G 
L(/i m+e (S 2 ), /i m+e (§ 2 )). Thus, by a standard perturbation result for infinitesimal generators of 
continuous analytic semigroups (see PQ and [M] ) , the desired assertion follows if we prove that 
-Qi(0) G H(h m+9 {§ 2 ),h m - 1+e {§ 2 )). 

Since (cf. [£]) 

= 1 - e[ri(u) + ^77(0;)] + o(e), 

where A w is the Laplace-Beltrami operator on the sphere § 2 , we have /Co(0) = /C(0) = 1 and 
£1(0)77 = — hAufj. Hence, from the definition of Qi we see that 

Si (0)77 = -v| S 2 • V^t/ - 7u^| g2 • v, (2.45) 

where v is the solution of the following boundary value problem: 

Act = f(o~) in \x\ < 1, 

V • v = g(a) in \x\ < 1, 

- Av + Vp - -V(V • v) = in |z| < 1, 
3 

a = 1 on \x\ = 1, 
T(v,p)n = —711 on \x\ = 1, 

v dx = 0, 

k|<1 

v x xci x = 0, 

'N<i 

and = — ^R(0)(A W 77 • n). Clearly, v = v s — the radially symmetric stationary solution of 
the problem (1.1)-(1.8) (see Appendix A), so that v|§2 = 0. It follows that 

Qi (0)77 = -7U„| S 2 • u = -ju v y • n = Jn • R(0)[n • A^t/] = -^Atf. (2.46) 



Define Aq = <9 n (A,tr g 2) 1 (0, •), i.e., Aqij = d n ip v for 77 G /i m+e (§ 2 ), where ip v is the solution of 
the following boundary value problem: 



= in |x| < 1, ip v = r] on |x| = 1. 

It is well-known that (cf. [S]), 
We rewrite 



(2.47) 



(2.48) 



Q 1 (0) = -1a -^(2A 1 -A ). 



In what follows we prove that 



2A l -A G L(h m+e (S 2 ),h m+e (S 2 )). 



(2.49) 



Note that if this assertion is proved, then the desired assertion follows. 

Since AqT] = dnip^ and A\j] = 2u^|§2 • n, we have (2A\ — Aq)t] = 4u^|§2 • n — <9 n Vv Since 
= — iR(0)(A w 7y-n), by definition of the operator R(0) we see that there exist q v G h m ~ 1+e (S 2 ) 
and (rj G M 6 such that 

£(0)(u v ,q v ,(v) = (0,0,-^^-11,0,0). 
Besides, a simple computation shows that 



£(0)(VVvO,0) = (o,0,T(VVvO)n, / V^dx, x xdx 

^ J\x\<l J\x\<l 



\x\<l J\x\<l 

0,0,2d n V^, / r)-ndS x ,Q), 



x=l 



Since = Aip v \ r=1 



1 d 2 ip v 2d^ v 1 



+ — jP + "o 



£^f^»2 ^^'Y* i^*2 



d 2 ^ 



r=l 



<9r 2 



r=l 



+ 29 n '0 r; + A^?7, we have 



9 



_ d 2 ^ 

r=l 9r 2 



, V7 /"^ 



r=l 



-A w r? • n - 2d n ip v ■ n + V ^(d^). 



Hence 



£(0)(4u„- V^,4g„,4C„) = (o, 0, 4d n ^ ■ n - 2V a ,(d n ^ ? ), ^ ^-ndS^O 
It follows that 

(2Ai - A )r? =4u r? |§ 2 • n - <9 n ^ = (4u^ - V^)|§ 2 • n 

=tr s2 { l 7£(0)- 1 (0, 0, 48 n ^ ■ n - 2V a ,(d n ^ ? ), y ^ r? • nd5 Xl o) } • n 

=-2R(0)(V a; (a n ^))| s2 -n + tr S 2{j r /:(0)- 1 (0,0,49 n ^-n,y' ^ • n^, 0) } • n 
=5i?y + B r]. 



It can be easily seen that Bq G L(/i m+e (8 2 ), h m+e (S 2 )). Furthermore, minor changes to the proof 
of Lemma A. 2 in [9J show that also B x G L(h m+e (S 2 ), h m+0 (S 2 )) (see Lemma B.l and Corollary 
B.2 in Appendix B for details). Hence (2.49) follows. This completes the proof of Lemma 2.6. 
□ 

Since H(h m+e (§ 2 ), /i m " 1+e (§ 2 )) is open in L(h m+e (S 2 ), h m - 1+e (S 2 )), from Lemma 2.6 we 
immediately get 

Corollary 2.7 For sufficiently small 6 we have 

- DQ(p) G H{h m+e (§ 2 ), h m ~ 1+e (§ 2 )) for p G Of +d {S 2 ). (2.50) 

By this corollary we see that, at least in a small neighborhood of the origin, the differential 
equation (2.39) is of the parabolic type in the sense of Amann pQ and Lunardi [24], so that 
the geometric theory for parabolic differential equations in Banach spaces presented in these 
literatures can be applied to (2.39). In the following sections we shall use this theory to prove 
Theorem 1.1. 

3 Linearization 

In this section we compute the spectrum of the operator DQ(0). Note that since h m+e (E> 2 ) is 
compactly embedded into /i m_1+e '(S 2 ), by Lemma 2.6 we see that the spectrum of the operator 
DQ(0) consists of all eigenvalues. 

To compute the eigenvalues of DQ(0) we first derive a useful expression of this operator. 
Consider a perturbation of the radially symmetric stationary solution (cr s , v s , p s , 0) (see (1.12)): 

a(x,t) = a s (r) + e<j)(r,uj,t), v(x,t) = v s (x) +ev(r,uj,t), p(x,t) = p s (r) + eij)(r,u,t), 

0(f) = {x G M 3 : r < l + eq(w,t)} (r = uj G S 2 ), 

where e is a small parameter, and <p, v( = (v%, V2, ^3)) , ip an d t\ are new unknown functions. 
From [18] and [19] we see that the linearizations of equations (1.1)— (1.8) are respectively as 
follows: 

A0 = /V-)0 in B 3 , (3.1) 
V -v = g'(a s )(f> in B 3 , (3.2) 

- Av + Vi;- -V(V-u) = in! 3 , (3.3) 
o 

(j) = -a' s (l)rj on S 2 , (3.4) 

T({7,V)n = -2 5 (l)V w ?? + 7(7? + -A u r))n + 4g(l)r/n on S 2 , (3.5) 

d t rj = v\ s2 ■ n + g(l)rj on S 2 , (3.6) 

/ vdx = 0, (3.7) 

J\x\<l 



vxxdx = 0, (3.8) 

x|<l 

Similarly as before, the system (3.1)-(3.8) can be reduced into a scalar equation in the 
unknown function 77 only. Indeed, given r\ £ C([0, oo), /i m+6l (§> 2 )), we first solve the second-order 
elliptic equation (3.1) subject to the boundary condition (3.4) to get <p(-,t) G /i m+e, (]B ) as a 
functional of 77, and next substitute this solution cf) into (3.2). It can be easily checked that 
(2.27) and (2.28) are satisfied by equations (3.2), (3.3), (3.5), (3.7) and (3.8). Thus by using 
Lemma 2.3 we get a unique solution (v(-,t),ip(-,t)) G h m ~ 1+e (M 3 ) x /i m ~ 2+6, (B 3 ) as a functional 
of rj. Substituting v = v{r,uj,t) obtained in this way into (3.6) and denoting 

B^T] = v\ s2 -n + 5(1)77, (3.9) 

we see that the system of equations (3.1)-(3.8) reduces into the scalar equation 

d tV = B lV . (3.10) 

Now, since the problem (1.1)— (1.8) is equivalent to the equation (2.39) with Q(p) given by (2.38), 
its linearization should correspondingly be equivalent to the linearization of (2.39) which reads 
as follows: 

dtV = DQ(0)rj. (3.11) 
Comparing (3.10) with (3.11), we get the following result: 

Lemma 3.1 DQ(0) = B~ ( . □ 

In the sequel we deduce the expression of in terms of Fourier expansions of functions on 
the sphere § 2 . 

For each I £ N U {0}, let Yi m {uj) (m = —I, — I + 1, • • • ,1 — 1,1) be a normalized orthogonal 
basis of the space of all spherical harmonics of degree I. Then {Yi m (u) : I = 0, 1, 2, • • • ; m = 
—I, — I + 1, ■ • • , I — 1, 1} is a normalized orthogonal basis of the scalar L -space on S 2 . As in [23], 
let Vi m {uj), Xi m (u>) and Wi m {oj), where / = 0, 1, 2, • • • and m = —I, — I + 1, ■ • • ,1 — 1,1, be the 
corresponding vector spherical harmonics. From [23j we know that all these vector spherical 
harmonics form a normalized orthogonal basis of the vector L 2 -space on § 2 (see also Appendix 
A of [18] and Section 2 of [22] for this assertion). Besides, for every IgNU {0} we denote 

d 2 2d l 2 + l 
1 dr 2 r dr r 2 

For simplicity of the notation, we shall not write out the whole expansions of (j), v, ip and 77, but 
instead merely consider each monomials in the expansions of these functions. This is reasonable 
because of the special forms of the operators appearing in the (3.1)-(3.8). Thus we put 

V = Yi m {uj). 

Then it can be easily verified that the corresponding solution of (3.1) and (3.4) is as follows: 



4>{r,u) = Fi(r)Yi m (u), 



(3.12) 



where Fi (r) is the unique solution of the following problem: 

L l F l (r) = f'(a s (r))F l (r) for < r < 1, F/(0) = 0, F,(l) = -^(1). 
Observe that by (3.2) and (3.3) we have 

A(V> - y((Ts)<l>) = in B 3 , 

Av - V( 5 '(<r s )0) = V(V> - p(<r s )(f>) in B 3 . 
Thus the solution of (3.2), (3.3), (3.5), (3.7) and (3.8) has the following expressions: 



(3.13) 

(3.14) 
(3.15) 



V>( r > w ) = ^9'(^s(r))4>(r,uj) + Pi m (r)Yi m (u) 

= ^Mr))Fi(r)Yi m (u) + P lm (r)Y lm (u) 



(3.16) 



V = a+b X X + Vi m Vi m + X lm Xi m + W lm W t 



Imi 



(3.17) 



where a, b are unknown constant vector, and Pi m (r), vi m (r), x/ m (r) and wi m (r) are unknown 
functions defined on [0,1] such that vi m {0) = xi m (0) = wi m (0) = 0. Using some well-known 
formulas for vector spherical harmonics (see [23] or [18], [22]), we have 



Vcj)-- 



F/(r) + ^(r) 



^? to ( w )+[iftr) + ^F,(f 



I 



21 + 1 



W lm (co), (3.18) 



+ 



^WMtM(t)] + ^'(^(r))^^) 



1 + 1 t? r ^ 



21 + 1 



2Z + 1 



rJ(r)tM + ^(r)# im ( W ), 



(3.19) 



Ay 

V(V-y5'(^)0) 



2Z + 



/ + \t-Wr Plm + l Plm] ^ lm+ \2T+T ^r Plm + (3 ' 20) 



V • v=- [v' lm {r) + lj ^ L vi m {r)] \j l^-Yi m (oj) + [w' lm (r) - l —^wi m (r)} J-^VT^ h,, U' )■ 



21 + 1 



(3.21) 



and 

Av = Li + i(vi m (r))Vi m (uj) + Li(xi m {r))Xi m {uj) + Li^ 1 (wi m (r))Wi m (uj). (3.22) 
By (3.15), (3.18), (3.20) and (3.22), we have 



i + i r 

21 + 1 



P[ m (r) + -P lm (r) = L l+1 (v lm (r)) - F}{r), 
r J 

Li{xim{r)) = 0. 



(3.23) 
(3.24) 



21 + 1 



By (3.14) we have 

By (3.2) and (3.21) we have 



PL(r) + —Plrn{r) = L^iwUr)) - F?(r). 



Li{Pi m (r)) = 0. 



I 



21 + 1 



I - 1 



Wlr. 



-wim{r) 



Tff^j v\ m {r) + l -^vi m (r) - . ( /(a,(r)).f',(/'). 



Solving the ODE problem (3.23)-(3.27), we get 

P lm (r) = 2(21 + 3)A ir l , x lm (r) = B x r l 



(3.25) 
(3.26) 

(3.27) 

(3.28) 



vi m (r) 



21 + 1 I + 1 (21 + 3) V 21 + 1 

r+1 2Z l+1 



' + 1 c/Wrmir) - r —— I s l+3 F}(s)ds 



-1-2 rr 

21 + 3 J 



21 + 11 + 1 



(3.29) 



and 



«•/,„</■) Civ 1 ' 1 + \l ^^-(2/ + 3)A ir l+1 



,. i i r l—l fRs 

- 5 '(<7 fl (r))FKr)- / s -^ F f( s )ds 



(21 - 1) V 21 + 1 



2/ - 1 



EE (V- 1 + 



^Ll(2l + 3)A 1 r l+1 -w l (r), 



(3.30) 



where Ax, B\ and Cj are constants. 

Next we consider the boundary condition (3.5). Again by using some well-known properties 
of vector spherical harmonics (see |23j or [18], [22]), we can rewrite v in (3.17) as follows 



v(r,u) = a + b x x + H a (r)Yi m (uj) u + H l2 (r)V w Yi m (u) 



(3.31) 



where 



21 + 1 



it • • I 1 + 1 f \ _l / 1 ( \ tj ( \ v lm(r) w im (r) 



21 + 1 



V(Z + 1)(2/ + 1) Vpl+l) 



Thus 



T(^,V)n= - 5 , (l)^(l)^ m (a;)^+ [2 J ff/ 1 (l)^ m (^)-V(l,^)]w 
+ - fli 2 (l) + H' l2 (l)]V u Y lm (uj). 

Note that by (3.16) and (3.28) we have 

V(l, u) = Y lm (u) [ - ^'(1)^(1) + 2(21 + 3) A,] . 



(3.32) 



Substituting this expression into (3.32) we get 



T(v, V)n = |</(1K(1) + 2H' n (l) + ~</(lK(l) - 2(2/ + 3)A 1 ]Y lm (u) u 
+ [H n (l) - H l2 (l) + H' a {l)]V u Yi m (u). 



(3.33) 



On the other hand, putting rj = Yi m (uj) in (3.5) and using the well-known relation A^Y^^u) 
-1(1 + l)Y lm (u), we get 



T(£,^)n = -2 5 (l)V w y m (cu) + [ 7 (1 



1(1 + 1), 



+ 4g(l)] Y lm (u)u. 



(3.34) 



Since V aJ Y5 m (u;) and Yi m (uS) lo are mutually orthogonal, by comparing their coefficients in (3.33) 
and (3.34) and using the relations 



iZll(l) = 

^i(l) = 
we obtain 



1 



21 + 1 



vi m (l) + 



I 



21 + 1 



w lm (l), H l2 (l) 



Vlm(l) 



Wl m (l) 

V(/ + i)(2i + i) ' #+!)' 



+ 



2T+\ vUl) + V 2TTT^ (1) ' 



<>(!) 

v/(i + l)(2/ + l) ' v/P+1)' 



+ 



\/^TT^ m(1) + V^TI^ m(1) = I (2 " /2 " l) + 2a{l) " 9 ' {l)<J ' s{l) + {21 + 3)Al ' 



and 



1 



1 + 2 m 

V2T+Ti-vm vlmil) + ^7f wlm{1) 



+ 



+ 



v /(/+i)(2/ + i) yipTi) 

We now proceed to consider the equation (3.6). By (3.9) and (3.29)-(3.31) we have 
B y Y lm (u) = v | g2 • n + 0(l)y Jm (w) = a • w + [iT,i(l) + s(l)]y, m (w) 
=a- w + y m (w) 5(1) - 



(3.35) 



2g(l). (3.36) 



1+1 
21 + 1 



Vlm(l) + 



I 



21 + 1 



=a • w + y m (cu) 5 (1) + l(Ar + C{) + 



2Z + 1 w V 2/ + 1 
where Ci = C\/ y/l(2l + 1). Thanks to the constrain condition (3.7) we see that 



, (3-37) 



_3_ 

47T 



/ {Hn(r)Y lm (io) u + H l2 {r)V u Y lm {w)} dx = for Z G N, I + 1 



(3.38) 



(cf. (5.8) in [18]). To compute A 1 + C x we substitute (3.29) and (3.30) into (3.35) and (3.36), 
which gives 



21 + 1 



1+-1^(1)=-1(2 -I 2 -I)- 2g(l) + g'(l)a' s (l) 
+Ct(l 2 - I) + (I 2 - I - 3)A U 



(3.39) 



and 

2/ 2 -I- 4/ 

= 2s(l) + 2d(l - 1) + -J^-M- (3.40) 
(3.39) x 2(2/ + 1) + (3.40) x 3(Z + 1) yields: 

Al + 6l = 2(/-l)(2P + 4 i + 3) {I (1 - ' )(2 ' 2 + 5i + 21 - < 4i + 2)9 ' (1 ^ (1 > 

+(2 '- 2)9(1)+ 7i|mf"'' i(1) -^ (4, - 1)i;(1) 



By the definitions of F\ [i = 1,2), v\ and u)j respectively in (3.20), (3.29) and (3.30), and by 
straightforward calculation we easily have 



= \l^l£^Mr))F,(ry + 'ir, (3.42) 



~ - t/aTT^ (1) ' (3 ' 43) 



= -^f^s'(iy,(l) - + f 9'Mrm(r)r>+* dr, (3.44) 

« 4 > = \/^TT27^T 9 ' (1) "- (1) - <3 - 45) 

Substituting (3.38) and (3.41) into (3.37) and using (3.42)-(3.45), we see that, for / / 1, 
g 7 y ^H= 2Z2+ 1 4Z + 3 {g(l)(2f + 3)(Z + 1) - Jf(2Z + 1)(Z + 2) 

+(2/ + 3)(Z + 1) ^ «/(<7 s (r)) J Fi(r)r ,+2 dr}y, m (u;). (3.46) 

We define, for / > 2, 

* = l(?+2K2l + lj I"* 1 ' + 1 ' ^.M W-V« *] , (3.47) 

Then in case ' > 2 (3.46) can be rewritten as follows: 

BjYi m (uj) = ai{i)Y lm {u). 
In the case I = we have, directly from (3.46), that (3.49) also holds with 

«0 = Ml) = 5(1) + f g'(cr s (r))F (r)r 2 dr. (3.49) 
Jo 



Finally, we consider the case I = 1. Since the problem (1.1)-(1.8) is translation invariant, 
by some similar argument as those in [6] and [18] we see that B^ij = for any sphere harmonics 
of degree 1. In particular, we have 

B ry Y lm (uj) = 0, m = -1,0,1. 
In summary, we have proved the following result: 

Lemma 3.2 DQ(0) = B~ is a Fourier multiplication operator having the following ex- 
pression: For any r\ G C°°(§ 2 ) with Fourier expansion r)(ui) = X]^o EL=-i c i™^m( w )» we ^ ave 

oo I 

B 7 7](uj) = aocooloo + ^ ^2 a i(Tf) c imYi m (uj), (3.50) 

1=2 m=-l 

where ao and o-iip/) defined in (3.49) and (3.48), respectively. □ 

As usual, for a given closed linear operator B in a Banach space X, we denote by p(B) 
and a(B) the resolvent set and the spectrum of B, respectively. The set of all eigenvalues of 
B is denoted by o~ p (B). As mentioned in the beginning of this section, we have o~(DQ(0)) = 
a p{DQ(0)). Hence, from Lemma 3.2 we immediately obtain the following result: 

Lemma 3.3 The spectrum of DQ(0) = B-y is given by 

o-iBJ = {a , 0} U {a/( 7 ) : I = 2, 3, 4, • • • }. 

Moreover, the multiplicity of the eigenvalue is 3. □ 

The next result shows some useful properties of olq and 7^ (I > 2): 

Lemma 3.4 We have the following assertions: 

(i) ao < 0. 

(ii) 7; > for all I > 2, and lim/^oo 7; = 0. 

(Hi) There exists an integer I* > 2 such that < 7/ for all I > I*. □ 

Proof: By Assumption (Al) we have /' > 0. Thus, by the maximum principle we see that 
Fo(r) < 0. Furthermore, since u(r) = —a' s (r) is a solution of the problem 

u"(r) + -u'{r) = f'{a s {r))u(r) + -^u for < r < 1, u(0) = 0, u(l) = ~cr' s (l), (3.51) 
by comparison we easily get Fo(r) < —a' s (r). Thus, since g' > (by Assumption (^42)), we have 

g'(a s (r))F (r)r 2 dr < - g (a s [r))a' s (r)r 2 dr = -g{l) + 2 / g(a s (r))r dr. 
Jo Jo Jo 

Hence, 

a = g(l)+ g'(a s (r))F (r)r 2 dr <2 g(a s (r))r dr < 0. 
Jo Jo 

The last inequality follows from the facts that g' > and / g(a s (r))r 2 dr = (by (A. 3) and 

Jo 

(A. 4) in Appendix A). This proves (i). 



Next, from (3.13) we have, for any I > 2, that 

F{'{r) + -F{{r)-- 9 F l (r)-f{o s {r))F l ir) 



I 2 + 1-2 



Fi(r) < 



Since /' > and u(r) = —a' s (r) satisfies (3.51), by comparison we get —a' s (r) < Fi(r) < 0. 
Hence 

g(l) + ! g'{a 8 {r))F l {r)r l+2 dr> g(l) - f g 1 ' (a s (r))a' s (r)r l+2 dr 
Jo Jo 

>g(l)- [ g'(a s (r))a' s (r)r 3 dr 
Jo 



g(a s (r))r 2 dr = 0, 



(3.52) 



so that 7; > for all I > 2. Moreover, since l-F^r)! < cr' s (l) and < g'(a s (r)) < <?'(!), we have 



ll 



4(2Z + 3)(Z + 1) 



g(l) + / g'ia^F^ry^dr] < 8Z" 1 [<?(!) + 



l{l + 2)(2l + l) 
Hence hmi_ >00 7; = 0. This proves (ii). 

Finally, by direct computation we have 

7m -11 = -%(!)(! + o{l))l- 2 as I oo. 
From this fact the assertion (Hi) immediately follows. The proof is complete. 
By virtue of the the assertion (ii) of the above lemma, we introduce 

7* = max 7;. 

Z>2 



1 + 3 



□ 



(3.53) 



Note that Lemma 3.4 ensures that < 7* < 00. 

It is interesting to compare the threshold number 7* defined above with the corresponding 
threshold number for the porous medium structured tumor model obtained by Cui and Escher 
[8], which we denote by 7*. Recall that 7* = max;>2 7;, where, for the case R s = 1 and a = 1, 



ll 



g(l)+ I g'^F^ry^dr 



for I > 2, I G N. 



(3.54) 



l(l-l)(l + 2 ) L J 

From Lemma 3.2 of [8] we know that {7;}z>2 has the same properties as {7/};>2 presented in 
Lemma 3.4. 

Lemma 3.5 7; < 7; for all I > 2, I G N ; so that 7* < 7*. 
Proof: For any I > 2 we have, by (3.47), (3.52) and (3.54), that 



ll ~ ll 



4(2Z + 3)(Z + 1) 



/(/-!)(/ + 2) l(l + 2)(2l + l) 



g(l)+ / ^{l)Fi(r)r l+2 dr 



< 



1(1 + 2) 
2 



1 2(2Z + 3)(Z + 1) 



1(1 + 2) 
This completes the proof. 



I 

1-1) 

□ 



21 + 1 



5 (1) + / <7'(l)F z (ry +2 dr 



g(l)+ / 5 '(l)^(ry+ 2 dr 



0. 



4 The proof of Theorem 1.1 



In this section we give the proof of Theorem 1.1. Note that since G a(DQ(0)), the standard 
linearized stability theorem for parabolic differential equations in Banach spaces cannot be 
applied to treat (2.39), and we have to employ the method of center manifold analysis. We shall 
construct a locally invariant center manifold, which consists only of equilibria, and show that 
this manifold attracts nearby transient solutions at an exponential rate. Similar method was 
applied in [7J, [E] and [26] . 

Proof of Theorem 1.2: We fulfill this proof in four steps. 

(i) By the definition of 7* we see that for any 7 > 7*, 0/(7) < for all I > 2. Besides, by 
(3.50) we see that is an eigenvalue of geometric multiplicity 3, and the kernel of DQ(0) = i3 7 
is the space X c := span{Y" lm ;m = —1,0,1}. Let X~- be the orthogonal complement of X c in 
L 2 (S 2 ), and for fixed m > 3 we denote hf +e (8 2 ) = h m+e (S 2 ) n X^. Then we have 

fc^(s J ) = fc^(s 2 )eJfc. 

This decomposition induces two projection operators 7r c and n s which map h m+s (E> 2 ) onto X c 
and /i™ +6, (§ 2 ), respectively. From Lemma 3.2 we know that /3 7 commutes with them. 

(ii) Let M(rj) = Q(rj) — £> 7 7/. Then the equation (2.39) can be rewritten as follows: 

d t rj = B^7] + M(rf) for t > 0, and 7/(0) = rj . (4.1) 

The little Holder spaces have the following well-known interpolation property 

(h ao (§ 2 ),h ai (§ 2 )), & = ft(i-*)°»+*n(S 2 ), if (1 - 0)a Q + 0<ti i Z 

where < < 1 and (•,•)■# denotes the continuous interpolation of Da Prato and Grisvard 
(see |24|). By Lemma 3.1 we know that i3 7 = DQ(0) generates a strongly continuous analytic 
semigroup on 

h m-l+9(g2j with domain h m+e (S 2 ). Thus by Propositions 6.2, 6.4 and Theorem 
6.5 in |10j we conclude that there exists an open neighborhood O of the origin in X c and a 
mapping 

c g c m (o, k™ +6 (§ 2 )) with c(o) = 0, ac(o) = 0, 

such that the 3-dimensional submanifold M c := graph (C) of h m+e (S 2 ) is a locally invariant and 
stable manifold for the evolution equation (4.1). Note that A4 C consists only of radial equilibria, 
i.e. Ai c is the set of all spheres of radius 1 with centers sufficiently close to 0. Furthermore, by 
the above-mentioned results of [10] we know that A4 C attracts at an exponential rate all small 
global solutions of (4.1) in /i m+6l (S 2 ). More precisely, there exists e > such that the solution 
to (4.1) exists globally for any 770 with ||?7o|L m + e (S 2 ) — e ' ana ^' moreover, there exist c > 0, K > 
and a unique zq = zq(t]q) G O such that for any t > there holds 

\\(ir c V {t),Tr s r){t)) - (^o,C(2o))||/ 1 -+«( S 2 ) < ^exp(-ct)||7r s 7 ?0 - C{7r c Vo )\\ hm+ e {s2) . (4.2) 

(hi) Now let 770 G h m+e (S 2 ) be given and ||?7o||c m + e (s 2 ) ^ £ - Then the solution of the 
equation (4.1) 77 G C([0, 00), h m+e (S 2 )) n ^((0, 00), h m+e ~ l (S 2 )), and it satisfies (4.1). By 
Lemma 2.1 and Lemma 2.5, it follows that the problem (1.1)-(1.9) has a global-in-time solution 



(a(;t),v(-,t),p(-,t),tt(t)), where Sl(t) = {x G R 3 : x = ru, < r < 1 + 77(0;,*), w G S 2 }. Since 
7W C is the set of equilibrium solutions which are sufficiently close to § 2 , there exists a xq G M 3 
such that (zo,C(zo)) = ^[x ]> where 77[ Xo ] is the distance function on § 2 introduced in Section 1. 
Then (4.2) implies that 

\\v(->t) - % ]IU m+e (s 2 ) ^ Kexp(-ct) for any t > 0. (4.3) 
By Lemmas 2.2-2.4 we have 

a(-,t) = K( V (t)) o <&-i } , v(-,t) = v(77(t)) o $-l )a p (.,t) = p(r/(t)) o c^, (4.4) 
Recalling the definition of (<T[ Xo ],V[ Xo ],P[ Xo ],tt[x ]), w e have 

<r[xo) = K(.V[ X0 ] ) %] 0] > U N] = v(T? [a;o ] ) o $-* o] , p [!Bo] = p(77 [a:o] ) o , (4.5) 

The explicit construction of (2.21) and the mean value theorem immediately imply that 
there is a positive constant C such that for any t > 0, 

||ft(77(t)) -ft(% ]) Hemes') ^ C 'll ? 7(*)-%o]ll cm +9(i 3 )- 
Then using (4.4), (4.5) we have 

\W(',t) -°"[xo])llc-+*(n(t)) ^ C||r/(*) -% ]||c-+e(n(t)) 5 ( 4 -6) 
for any t > 0. Similarly by Lemma 2.3 and Lemma 2.4 we also have 

||v(-,i) - V[ X0 ])|| C m-i + 9 ( Q W) < C||r/(t) - r?[ a;o ]|| c .m-i+e(n( t )), (4.7) 

l|p(->*)-P[ao])llc7m-3+fl(fl(t)) < CIW*) ~ %o] ||cf*»-»+«(ft(*)) > ( 4 ' 8 ) 

for any t > 0. Combining (4.2), (4.3) and (4.6)-(4.8), we see that (1.14) holds. 

(iv) Finally, if < 7 < 7* then by Lemma 3.3 we see that c(£> 7 ) H {A G C : ReA > 0} is not 
empty. It follows from Theorem 9.1.3 in [23] that the zero equilibrium of (4.1) is unstable. This 
completes the proof of Theorem 1.1. □ 

Appendix A: Radially symmetric stationary solution 

Theorem A Under Assumptions (Al)-(A3), the problem, (1.1)— (1.8) has a unique radially 
symmetric stationary solution (a s ,v s ,p s ,£l s ) with components having expressions in (1.12). 



Proof: It suffices to consider the following problem: 

<W + -°' s {r) = f(a s (r)) for < r < R„ (A.l) 
r 

a' s (0) = 0, a s (R s ) = l, (A.2) 

v' s {r) + -v s {r)=g{a s {r)) for < r < R s , (A3) 
r 

v,(0) = 0, v s (R s ) = 0. (A4) 



[|a + [Psir) ~ \9Mt))] =0 for < r < R s , (A5) 

K(0) = 0, p a (R a ) = ^- + p(l). (A6) 

Indeed, it is straightforward to verify that for a radially symmetric stationary solution the 
equations (1.1) and (1.2) have respectively the forms (Al) and (A3), the boundary condition 
(1.4) has the form of the second equation in (A2), and (1.6) becomes the second equation in 
(A4). Next, taking divergence in both sides of (1.3) and using (1.2), we see that (A5) holds, 
by which we have 

T(v s ,p s )n\ r=Rs = [2v' s (R s ) - P s(Rs) ~ |s(l)]n 
= [2g(l)-p s (R s )-^g(l)]n 
= [^9(1) - Ps(Rs)]n, (by (A3), (A4)) 

where n(x) = x/\x\ for x £ M 3 \{0}. Hence, the second equation in (A6) follows from the 
boundary condition (1.5). Finally, the first equations in (A2), (A4) and (A6) are imposed 
to rule out possible solutions possessing singularities at r = for the problem without these 
equations, which are not meaningful solutions of (1.1)-(1.8). 

From [4] we know that under Assumptions (Al)— (A3), the problem (Al)-(A4) has a unique 
solution (a s (r),v s (r), R s ). Besides, using (A5) we immediately see that the function 

7 4 
Ps(r) = — + -g(v s (r)) 

solves (A5) and (A6). Since the solution of (Al)-(A6) is obviously unique, we see that the 
desired assertion follows. □ 



Appendix B: Boundedness of the operator B\ 

Lemma B.l Let m e N, m > 3, and < 9 < 1. Let h £ (/i m ~ 2+e '(S 2 )) 3 satisfy h • n = 
and let v = R(0)h. Then tr§2(v) • n £ /i m+e (§ 2 ), and there exists a constant C independent of 
h such that 

||tr S 2(v) • n|| ftm +e(§2) < C||h||( ftm - 2 +e(§2))3. 



Proof: Fix a function d £ BUC°°(R 3 ) such that 

tr §2 (d) = 0, tr g2 (Vd) = n, d n (Vd) = 0. (5.1) 

Extend the normal vector field on § 2 into M 3 by setting n = Vd. From the proof of Lemma 
2.3 we know that there exists a unique (v,p,C) £ (/i m_1+e (S 2 )) 3 x /i m ~ 2+e (S 2 ) x M 6 , where 
v = (y±, V2, V3), such that 

£(0)(v,p,C) = (0,0,h,0,0). (B.2) 



We have 

Ap = V • (Av - l (Q) = in B 3 . (B.3) 

Since, on one hand, 

T(v,p)n • n = 2niiijdiVj — p = 2njd n Vj — p 
= 2<9 n (v • n) - 2v d n (Vd) - p 
= 2<9 n (v ■ n) — p on § 2 , 

and, on the other hand, by (B.2), 

T(v,p)n • n = h • n = on § 2 . 

Hence, 

2<9 n (vn)-p = on S 2 . (B.4) 
Define * = v • n — \pd. Then we have tr§2(v) • n = tr§2(\I/) and 

A^ = A(v -n-\pd) 

= Av • n + 2d i v j d i n j + v • An - ±pAd -Vd-Vp- \dAp 
= (Vp + l (()) • n + 2d i v j d i n j + v • An - ±pAd - Vp • n 
= Z (C) • n + 2diVjdiYij + v • An - \pAd in B 3 . 

Furthermore, by (B.l) and (BA) we have 

= 3 n (v • n) - ^d n p - |n • (Vd) = 3 n (v ■ n) - X -p = on § 2 . 
Hence ^ is the solution of the problem 



A* = l (() ■ n + 2diVjdinj + v • An - \pAd in 
d n ^ = on S 2 . 



(S.5) 



It follows by classical Holder estimates for second-order partial differential equations of the 
elliptic type that 

||tr §2 (v) • n|| ft m+e( §2 ) = ||*L™+e(§2) < CH^II^m+^i 3 ) 

<C(||A*|| /im _ 2+e(S 3 ) + ||^*|| ftm -i+« (sa) + ll^H^-a+s^Sj) 

<C(||v|| ftm _ 1+e(S 3 ) + ||p|| /im _ 2+e(5 3 ) + ICI) 

^C||h||(/i m - 2 +"(s 2 )) 3 - 

The proof is complete. □ 

Since V u} (d n ip v ) • n = uj ■ V w (9 n V' ?? ) = 0, by Lemma B.l we immediately obtain 

Corollary B.2 Let B x r) = -2R(0)(V w (<9 n ^)) | g2 • n for rj £ h m+e (S 2 ). Then we have 
B\ G L(h m+0 (E> 2 ),h m+9 (§> 2 )). a 
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